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The influence of an inhomogeneous screened electric field on the charging of dust grains in a 
thermal dusty plasma is studied. The electric field of charged grains is considered within the cell 
approach, where the problem is reduced to a one-particle one. Within the model of quasichemical 
equilibrium, which is generalized to the case of an inhomogeneous screened electric field, we obtain 
the value of mean charge of dust grains, the distribution function of grains over charges, and the 
variance of this distribution. We also give the criterion of an inhomogeneity g[z) of the electric 
field and show that the infiuence of screening effects can be neglected in the case of the rarefied 



subsystem of dust grains (for r^/r^ ^ 1) and in the case of small-radius grains (rp 
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I. INTRODUCTION 

The thermal dusty plasma characterized by the equal- 
ity of the temperatures of all components occurs mainly 
under the terrestrial conditions in the combustion prod- 
ucts of various fuels One of the important problems 
of the theory of a dusty plasma is the problem of deter- 
mination of the charge of dust grains 0,13 • This problem 
is solved with the use of various methods considering its 
specificity with that or other degree of adequacy. Below, 
we discuss the most typical works. 

In work the existence of an analogy between the 
process of thermoionization of dust grains and the pro- 
cess of ionization of atoms was used for the determina- 
tion of the charge of dust grains. The charge of atoms at 
the ionization varies by at most several values of electron 
charge. Therefore, it is possible to neglect the screening 
effect and to consider that the electrons pass a sufficiently 
large (infinite) path, by leaving the ions emitting them. 
The same approach was applied in work Q to the deter- 
mination of the charge of identical dust grains of radius 
Tp. Within this method, the analog of the ionization 
energy / in the Saha equation is the work function of 
electrons leaving the surface of grains, which have charge 
z: I => Wo -|- ze^/rp, where Wq is the work function for 
the surface of neutral grains, and the term ze^ jr^ [z is 
the charge in units of the electron charge e) has meaning 
of the work needed to move an electron from the surface 
of a grain to infinity. It was shown that, under the made 
assumptions, the mean charge Z of dust grains for Z 3> 1 
is given by the expression 



kTrp Tics , 1 
Z = — TT- In 1- - , 



(1) 



where Ties = 2(27rmofcr/ft,2)3/2 

exp{—Wo/kT) is the equi- 
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librium density of electrons near the surface of an emit- 
ting grain, and ric is the mean density of electrons in a 
plasma. 

Within the approach developed in which will be 
called quasichemical, works [a-Q considered the possibil- 
ity to form a negative charge and constructed the distri- 
bution function over charges of grains. It was shown that 
the coefficient kTr^/e^ is the distribution variance. We 
note that a small number of dust grains can be negatively 
charged, but, on the average, the dust grains in a plasma 
of combustion products have a positive charge. 

In work Q , the charge of dust grains of the same size rp 
is determined from the condition of balance of currents: 
the thermoemission current is balanced by the current 
of electrons, which recombine with charged dust grains. 
The thermoemission current is determined, like that in 
U , by the Richardson-Dushman equation. It is assumed 
that all electrons colliding with a dust grain recombine. 
The recombination rate is estimated in the approxima- 
tion of binary collisions as — ancJip, where a = av, a is 
the cross-section of a dust grain [a k, Trrp), v is the mean 
thermal velocity of relative motion, and rip is the mean 
density of dust grains in a plasma. The balance equation 
for flows allows one to establish the connection between 
the mean charge Z and the mean density of electrons n^: 
Z = kTr^l Innes/j^e- In fact, this result coincides with 
([T]). In the same work, the value of Uc was taken from 
experimental data (in principle, it can be calculated theo- 
retically), which gave the estimate of the value of charge. 
For example, for ric = 3.47 x 10^° cm~^, rip = 2 x lO'^ 
cm~'^, T — 3285 K, and for the grains with the radius 
Tp = 10~^ cm, the charge Z was calculated to be 1700. 

There is no basic difference between the above- 
mentioned approaches, which was indicated in [Toj . 
where the formula of equilibrium constants was obtained 
by the kinetic way. This was further used in [3~Q|- 

Works [Tl| - {l5| present the attempts to consider the in- 
homogeneous distribution of electrons near a dust grain 
and the influence of the appropriate electric field on the 
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recombination. In [lT| , the screening effects were consid- 
ered in the hnear approximation. Work [l5| dealt with a 
separate dust grain, and the distribution of the electric 
field in its neighborhood was found as a solution of the 
Poisson equation in the Debye-Hiickel approximation. 

A certain contradiction is inherent in the above ap- 
proaches: the introduction of the notion "the mean den- 
sity" is not proper, because it was assumed that the ther- 
moemission electrons go to infinity. 

Therefore, some different approaches were naturally 
developed. There, it was considered that the thermoe- 
mission electrons remain in a bounded region including 
the dust grain [l6l - l20| . This region (we will call it the cell) 
is electrically neutral. The electrons can pass from one 
cell to another one, but the mean number of electrons in 
each cell under the given conditions remains fixed. Such 
an approach allows one to consider only one single dust 
grain in a cell, i.e., the problem becomes one-particle. 
We note that, in the frame of cell approaches, the intro- 
duction of the mean density of emitted electrons is quite 
proper. 

The system composed from identical spherical dust 
grains, electrons emitted by them, and a buffer gas, was 
considered in [31 . Each dust grain is located inside an 
electrically neutral spherically symmetric cell of radius 
Tc. The size of a cell is determined from the geometric 
reasoning as half the mean distance between dust grains: 

The dust grains have the same charge Ze. The electric 
field inside a cell is determined by the numerical solution 
of the Poisson equation under relevant boundary condi- 
tions. In (T6| . it was accepted that, on the boundary of 
a cell, the potential and the strength of the electric field 
are zero. For the cases of weak and strong screenings of a 
grain by the electron cloud, the approximation formulas 
for the electric field potential were obtained. In work [T7| . 
it was shown that the solution of the electrostatic prob- 
lem in the case of the weak screening (rD ^ ''p) transfers 
in ([1]), if the subsystem of dust grains is rarified. 

In (Tsl . ITgj . the electric field inside a cell was deter- 
mined by the solution of a linearized Poisson equation 
supplemented by two boundary conditions, which estab- 
lish: (i) the equality of the electric field strength on the 
boundary of a cell to zero and (ii) the connection of the 
electric field strength on the surface of a grain with its 
charge. 

We note that the boundary condition (i) is trivial (it 
follows from the Ostrogradskii-Gauss theorem) and gives 
no new results. Moreover, the influence of the electric 
field on the charging of dust grains was not considered in 
the proper way. 

Work [13, which was also based on the cell approach 
to the description of properties of a dusty plasma, used 
the notion "the plasma potential" ("bulk plasma poten- 
tial") introduced earlier in [2l|, [2^ . It was proposed to 



reckon the electric field potential from this "plasma po- 
tential". The authors commented on such a choice of the 
reference point for the potential that, only in this case, 
one can set 0(r)|r->.oo — ^ 0. In the previous works, "the 
plasma potential" was defined so that the Boltzmann dis- 
tribution was consistent with the condition of the electric 
neutrality of a plasma. This proposition is not proper, 
since the process of thermoelectronic emission from an 
isolated grain is nonstationary, and the use of the Boltz- 
mann distribution is impossible. 

The consideration of the processes of charging and the 
influence of the screening effects on them has become es- 
pecially actual after the discovery of ordered structures 
in the thermal dusty plasma |23| . This problem was ana- 
lyzed, in particular, in works |24| - |26| . where the influence 
of screening effects on the charging of a separate grain 
was studied. In work [23|, the behavior of the screened 
field of a dust grain, which was described with a linearized 
Poisson equation, was studied by numerical methods. It 
was shown that the essential deviation from a linear the- 
ory of screening should be expected for grains, whose 
radius is of the order of the Debye one. The dynamics of 
the charging of a dust grain in the presence of external 
sources of ionization with regard for the photoemission 
from its surface was numerically studied in [25| . The 
value of charge was determined from the condition of 
equality of flows. The densities of electrons and ions in 
the vicinity of a dust grain and the charge of a grain as 
a function of the time were calculated. In work [26|, the 
influence of boundaries on the screening of a point-like 
dust grain was studied. 

In addition, the Saha equation became again used as 
a tool to describe the properties of a dusty plasma in 
the recent works [13, [Hi (see also [U, [H]). This is ex- 
plained by the simplicity and the physical clearness of 
this approach. Therefore, the necessity to comprehen- 
sively study the possibility of the use of the quasichemical 
approach seems obvious. 

In the present work, we study the value and the vari- 
ance of the mean charge of dust grains in the thermal 
plasma. To this end, we use: (i) the cell approximation 
and (ii) an approach based on the quasichemical model of 
the charging of dust grains. No buffer gas is considered. 



II. DEFINITION OF A CELL AND THE 
ELECTROSTATIC FIELD ENERGY 

Here, we will formulate the basic positions of the cell 
model and will determine the energy of the electrostatic 
field created by a charged dust grain. The limiting case 
of a weakly charged grain will be analyzed as well. 



A. The cell model of a dusty plasma 

Let us consider the identical dust grains of radius rp 
with the mean charge Ze. which are in equilibrium with 
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electrons emitted by them. We assume that it is possible 
to separate an electrically neutral cell around each grain. 
The cell radius rc is given by relation ([2]) . The condition 
of electric neutrality of a cell takes the form 



Ze + ^TT / p{r)r^dr = 0, 



(3) 



where n^o is the mean density of electrons for = 0, 
i.e., on the boundary of a cell. We consider that the 
screening of a grain is realized only by electrons emitted 
from its surface. 

Equation Q is supplemented by boundary conditions 
on the boundary of a cell and on the surface of a dust 
grain: 



where the distance r is reckoned from the center of a 
grain. 

We assume that the distributions of the bulk charge 
p{r) and the potential (p{r) of the electrostatic field inside 
a cell have spherical symmetry: p{r) => p{r) and (p{r) 

The distribution of the potential is described in the 
approximation of self-consistent field: the potential 0(r) 
satisfies the Poisson equation, in which the charge density 
p{r) is determined by the Boltzmann distribution. 

We note that the charge density in the vicinity of a 
charged grain is not always described by the Boltzmann 
distribution. For example, it is not proper if the nonsta- 
tionary problems, involve a change of the charge of a dust 
grain in the course of time. In this case, it is expedient 
to use the methods of no neq uilibrium thermodynamics, 
as it was proposed in 0, l30l |. 

In order to solve the problem posed in Introduction, we 
turn to the linearized Poisson equation. This is justified 
by the following reasoning. The potentials, which are the 
solutions of the linearized and nonlinear Poisson equa- 
tions, are close for rp < r < rp + td and rp + td < r < Vc 
(rD is the Debye radius of dust grains). In the second 
region, the linearized Poisson equation is generally quite 
adequate. The considerable difference between the po- 
tentials is significant only for r ^ td. This circumstance 
is of importance in the problems, in which the detailed 
behavior of the potential is defining. To estimate the 
mean charge Z within the method used by us, it is nec- 
essary to calculate the energy of the electrostatic field, 
which is a functional of the potential. Therefore, the 
"fine" details of the behavior of the potential are insignif- 
icant. 

In the dimensionless form, the linearized Poisson equa- 
tion reads 



1 



X 



0, 



(4) 



where Af is the dimensionless radial part of the Laplace 
operator. The dimensionless parameters are as follows: 



e(/)(f ) 
kT ' 



(5) 



where c and 'ijj{f) are, respectively, the dimensionless ra- 
dius of a cell and the dimensionless potential of the elec- 
trostatic field inside a cell, e is the elementary charge, 
and k is the Boltzmann constant. 

The quantity A is the dimensionless Debye radius equal 

to 



A 




kT 



(6) 



dip{f 



df 



Z 



(7) 



where = kTrp/e^ (for grains of a radius of 10 ^ cm 
at the temperature T = 3000 K, Zo ^ 10^ by the order 
of magnitude). The solution of Eq. (|4]) satisfying the 
boundary conditions ([7]) takes the form 

1 {Z/Z,)\ sh + , (A sh + ch f^) 
' f Ash^ -f ch^ ' ^ ' 

It is possible to show that the electric field strength 

l(Z/Zo)chi^-c(ish^+ch£^) ^ 



Ashi-i 



chi^ 



1 (Z/Zo)Ash^+^(Ash^+ch^^) 



Ashi^ 



chi^ 



(9) 



(here, E{f) = eE{r)r-p/kT) on the boundary of an elec- 
trically neutral cell is equal to zero. This is a natural 
consequence of the Ostrogradskii-Gauss theorem. Hence, 
the chosen boundary conditions ([7|) are quite suitable for 
the complete solution of the electrostatic problem, and 
the boundary condition E{rc) = should not be consid- 
ered instead of (j){rc) = 0. We note that i?(^) = with 
regard for the dependence of the Debye radius A on the 
grain charge Z (see the following subsection). 

The energy of the electrostatic field Wc\ can be calcu- 
lated in the standard way: 



Sir J 

V 



(10) 



Performing all necessary calculations, we obtain 
that the dimensionless energy of the electrostatic field 
Wci{Z) = Wci{Z)/kT reads 



W,i{Z) = aZyZo-f3Z + ^Zo, 

where the coefficients a, /3, and 7 are as follows: 

Ash2^ + 2A2 ch2^ + 2(<r - A^ - 1) 

a = ^ -. 

8D 

^ 2XD 



(11) 
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Figure 1: Debye dimensionless radius X{Z) for various radii 
of a cell <; 
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Figure 2: Schematic distribution of the density of thermoe- 
mission electrons in an electrically neutral cell. Here, Ties is 
the equilibrium density of electrons near the surface of a grain 
{r — Tp), Tioo is the density of electrons for tf) = {r — Tc), 
and rio is the mean density of electrons in a cell 



Since n^o = A^ , Eq. (|14p gives also the dependence 
of the density rico of thermoemission electrons on the 
boundary of a cell on the mean charge of a grain Z . The 
knowledge of the value of nco allows us to determine the 
mean density rio of thermoemission electrons in a cell: 



/■ / ^ 2 , 
— / Tip r r dr, 



(15) 



D = ( Ash^-^ +ch'^ ^ 



where Vc = 47r(rp — ?'p)/3 is the cell volume free of a dust 



A 



A 



(12) grain. In view of nc{r) ~ Uco ^1 + ^^^^ 



we obtain 



These coefficients depend implicitly on the charge of a 
dust grain Z. To take this fact into account, we deter- 
mine the dependence X{Z). 



B. Dependence of the Debye radius on the charge 
of a dust grain 

It is easy to verify that the equation of electric neu- 
trality ^ can be written in the form 



Z 



1 



ip{f)f df = 0. 



(13) 



By substituting solution ^ in p^ . we obtain 



^ = i [(, - A^) sh + A(, - 1) ch i^] . (14) 

This equation establishes the interrelation between the 
dimensionless Debye radius A and the mean charge Z of 
a grain. The corresponding dependence for various sizes 
of a cell <; is presented in Fig. [1] It is seen that, as the 
size of a cell increases (i.e., as the mean distance between 
dust grains increases), the dimensionless Debye radius A 
becomes much more than 1. 



nc{Z) = 



Z 



IZn 



Air 



(16) 



It is clear that the calculated mean density of electrons 
Tie in a cell in the limiting case oo should coincide 
with TicQ (see Fig. (2]), which has meaning, in this case, of 
the mean density of electrons in the system. 

Equation (jl4p is significantly simplified in the limiting 
case of small charges of a dust grain {Z ^ 0), which cor- 
responds to the condition A — oo. In this case, we have 
a small parameter y/fico and can expand the right-hand 
side of Eq. (HH) in the series in it. To within (Z/Zq)^/^, 
we obtain 



fieoiZ) 



Z 



02^0 



1 _ "4 Z 

a? Zq 



(17) 



where 02 = (-^^ - l)/3 and 04 = (<?^ - 5<;^ + 5<?2 - l)/30 
are the coefficients of the relevant degrees of y/heo- The 
domain of applicability of the asymptotic formula ([T7| is 
determined by the inequality hco 02/04. 
Thus, as A — > 00, we have the relation 



Zo 
ao — 
Z 



1/2 



HZ) 

which is valid for A 3> \/ 02/04 



1 + ^ — 

2alZo 



(18) 
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C. Behavior of the potential and the energy of the 
electrostatic field for large Debye radii (A ^ 1) 

In this limiting case, the obtained formulas for the 
potential ^ and the energy ([TT|) of the electrostatic 
field are considerably simplified. We now show that, as 
X oo, we arrive at the case considered in Appendix. 

Indeed, the condition A — > oo is equivalent to Z — > 0. 
Therefore, taking ^5)1 into accoimt and expanding ([5]) in 
the series in Z, we obtain (jA3p . Such a behavior of the 
cell potential was indicated in work fv^ . 

If we restrict ourselves by the first term in the expan- 
sion, then the formulas for a, /3, and 7 become 
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Figure 3: Function g{z) for various values of the dimensionless 
radius of a cell ? 
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Hence, the electrostatic energy (jlip in the approximation 
A — > 00 (Z — > 0) is equal to (jA5|) . In the limiting case — > 
00, the energy Wd passes into the energy of a nonscreened 
grain: ^ {Zef/2rp. 



III. THE MODEL OF QUASICHEMICAL 
EQUILIBRIUM 



As was mentioned above (see Introduction), the factor 
^ig'^ ze'^/rp in relation (|2ip represents the additional work 
made by an electron, by moving from the surface of a 
grain with charge z to infinity. 

The consideration of the influence of the electric field 
becomes essential in a number of cases. For example, for 
grains of the radius Vp = 10~^ cm with Z = 100, the ratio 
of the factor Ze^/rp to the work function Wq = 2.75 eV 
is equal to 0.52. We obtain the same value also for grains 
of the radius rp = IQ-^ cm with Z = 1000. Thus, the 
influence of the electric fleld becomes essential (for grains 



In this section, we recall the main results of works 0- 
01 and generalize the approach developed in them. For 
this purpose, we use the results obtained in the previous 
section. We will construct the distribution function of 
grains over charges, like that in but with regard for the 
performed corrections, and will determined the variance 
Zd of this distribution. 



of the radius rp = lO^^-lO"* cm) for Z - 100-1000. 
Such charges are typical of a dusty plasma. 

The Saha equation rip^'no/np^ ^"^ = Kz allows one to 
represent the mean density of dust grains with charge z 
as rip^-* = W^iZa^z/nc and the mean densities of dust 
grains Up and thermoemission electrons n^. as 



E 



E 



(22) 



A. Generalization to the case of a screened field 

It is known |31j] that, from the thermodynamic view- 
point, the ionization equilibrium is a particular case 
of the chemical equilibrium corresponding to simultane- 
ously running "reactions of ionization." These reactions 
can be written as follows: 



Ai = A2 



(20) 



where the symbol Ao means a neutral atom, Ai, A2, ... 
are atoms ionized one, two, etc. times, and e~ is an 
electron. 

Under the assumption that such a model of charging 
holds also for dust grains, the authors of works 0^ ob- 
tained the equilibrium constants 



Tics CXp 



z e 

TpkT 



(21) 



In the above-mentioned works, the condition of electric 
neutrality rie = Zup yields the formula for the mean 
charge Z of dust grains: 

CJO 

E ^exp 



(- 



2r„kT 



;ln^ 



00 / 

E exp ( - 

— — rv-j ^ 



z(z-l)e^ 
2r„fcT 



;ln2£^ 



(23) 



Equation for Z is not closed, since the quantity 
Tio remains undetermined. In works 0-0], it was taken 
from experimental data. In addition, the probability for 
a grain to obtain a negative charge is determined, obvi- 
ously, only by the energy of the electrostatic interaction. 
Therefore, for z ^ 0, we should set Ties = nc and we 
obtain Innos/jT-c = in (j23p . 

We now generalize the result in by considering 

(i) the existence of a finite domain (electrically neutral 
cell); (ii) the infiuence of screening effects. Then the 
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equilibrium constants (pi]) take the form 

W,^{z)-W,^{z-l) 



if. 



Ties exp 



kT 



(24) 



where the difference Wci(2;) — Wc\{z — 1) represents the 
additional work made by an electron, by moving from the 
surface of a grain into the cloud of electrons surrounding 
a grain. 

In Fig. [21 we present the function 



Wel(z) - Wel(z - 1) . , 

27 = ff(^) 



(25) 



for various sizes of a cell (t. It is seen from the fig- 
ure that, for ^ 3> 1, i.e., for distances between dust 
grains to be much more than their radii, the function 
g{z) — > 1. Hence, we may consider that the thermoemis- 
sion electrons move from grains to infinity. Therefore, 
for Tc S> ^p, we may neglect the screening and use the 
equilibrium constants (PT|). 
The equilibrium constants 



lead to the result 
E zexp[/(z,Z)] 



E exp[/(z,Z)] 



(26) 



where 



-W,,[z), z ^ 0; 

-W'ei(^) + zln^, z>0. ^ ' 



Here, the mean density of thermoemission electrons 
in a cell is a function of the mean charge of dust grains 
Z (see Eq. p^). 

The mean charges of dust grains Cc02 (Wq = 2.75 eV) 
for various temperatures and sizes of a cell (or, what 
is the same, for various mean distances between dust 
grains) are given in Tables IB lHIl 

First, we calculated the charge Z for given c, T, and 
Tp by formula (|26p (column Zcou)- In this case, we de- 
termined self-consistently the values of mean density of 
electrons rie in a cell by formula Then we calcu- 

lated the values of charge by the "old"' formula (col- 
umn Zoo ) , where we used the calculated value of rie . The 
computations were performed with the help of the devel- 
oped algorithm in language Java with the use of class 
BigDecimal from packet Math. 

In Fig. [31 we present the mean charge Z of dust grains 
Ce02 of the radius rp = 5 x 10~^ cm versus the dimen- 
sionless radius of a cell at T = 2500 K. It is seen that 
the role of the inhoniogeneity is significant for small cells. 
As the size of a cell increases, formula passes in 
In addition, the value of Z decreases for ^ ^ 5 and starts 
to grow for > 5. We explain such a result by that the 
quasichemical approach considered in the present work is 
not applicable for very dense systems (^ ^ 1). 
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Figure 4: Mean charge Z of dust grains Ce02 of the radius 
r-p = 5 X 10^''' cm vs the dimensionless radius of a cell ? at 
T = 2500 K 



B. Distribution of grains over charges 

The undoubtful advantage of the approach under con- 
sideration is the possibility to construct the distribution 
function of dust grains over charges like that in work Q , 
where the equilibrium constants (|2ip were used. Thus, 
by assuming that the charge of grains varies continuously, 
the following relation was obtained: 



lip 



2Zn 



<3^P 1 TT^ — " + -^^^^ ) = Pli^)- 

(28) 

The variance of this distribution is Zq = kTr-p/e^. 

In the cell approximation, the distribution law for the 
charge z of dust grains reads 



it,p 



exp [f{z,Z)] 

OO 

E exp[/(z,Z)] 



P2iz). 



(29) 



In Fig. [51 we present the distribution functions pi{z) 
and p2{z) at T = 2500 K, rp 5 x 10"^ cm, and ^ = 15. 
It is seen that the insignificantly small part of grains is 
negatively charged. 

The variance Zd of distribution (j29p can be determined 
in the standard way as 



Zu = P2 — Ml, 



(30) 



where fii is the i-th moment of the distribution, which is 

OO 

fH= ^>2(^). (31) 



Z — ~QO 



For the same parameters as in Fig. [4l we give the depen- 
dence Zd . It is seen that the variance Zd is somewhat 
different from Zq for small radii of a cell, and ^d('^) ^ Zq 
as (T — >■ OO. 
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Table I: Mean charges of dust grains calculated by the generalized formula (|26[l and by formula (|23p for various radii of a cell 
<; and grains Vp at the temperature T = 2000 K 





T = 2000 K 


rp = 10 * cm 


rp = 5 X 10 ^ cm 


r-p = 10 ^ cm 


■^coll 


rieo, cm^^ 




ZccU 


noo, cm"'* 


Zoo 




TicO, cm"'* 




5 


757 


1.46 X 10^2 


427 


261 


4.02 X 1012 


153 


13 


2.55 X 10^^ 


9 


10 


863 


2.06 X 10" 


662 


338 


6.46 X 10" 


263 


30 


7.10 X 10^2 


24 


15 


942 


6.67 X 10^° 


797 


385 


2.18 X 10" 


328 


40 


2.84 X 10^2 


35 


25 


1061 


1.62 X 10^° 


966 


449 


5.49 X 10^° 


410 


54 


8.24 X 10" 


50 


35 


1149 


6.40 X 10^ 


1078 


495 


2.20 X 10^° 


465 


63 


3.53 X 10" 


60 


50 


1248 


2.38 X IQS 


1196 


545 


8.33 X 10" 


523 


74 


1.41 X 10" 


71 


150 


1582 


1.12 X 10« 


1562 


713 


4.04 X 10* 


705 


107 


7.57 X 10^ 


106 



Table II: Mean charges of dust grains calculated by the generalized formula (|26p and by formula (|23|) for various radii of a cell 
<; and grains rp at the temperature T = 2250 K 



? 


T = 2250 K 


r'p = 10 * cm 


rp = 5 X 10~^ cm 


rp = 10 ^ cm 




rieo, cm"^ 




ZccU 


rioo, cm~^ 


Zoo 




rioO, cm~^ 




5 


1250 


2.41 X 10^2 


676 


472 


7.27 X 10i2 


264 


36 


7.01 X 10^^ 


23 


10 


1265 


3.02 X 10" 


955 


520 


9.95 X 10" 


398 


56 


1.35 X 10^2 


45 


15 


1326 


9.38 X 10^° 


1113 


562 


3.18 X 10" 


474 


68 


4.81 X 10^2 


59 


25 


1441 


2.20 X 10^'' 


1308 


627 


7.66 X 10^° 


570 


83 


1.27 X 10" 


77 


35 


1533 


8.54 X 10^ 


1436 


675 


3.01 X IQio 


633 


94 


5.28 X 10" 


89 


50 


1641 


3.14 X 10'' 


1571 


731 


1.12 X IQio 


700 


105 


2.01 X 10" 


102 


150 


2013 


1.42 X 10* 


1988 


918 


5.20 X 10* 


907 


143 


1.01 X lO^'' 


142 



Table III: Mean charges of dust grains calculated by the generalized formula (|26p and by formula (|23p for various radii of a 
cell <; and grains Vp at the temperature T = 2500 K 



f 


T = 2500 K 


Tp = 10 cm 


rp = 5 X 10 ^ cm 


rp = 10 ^ cm 




rieO, cm"^ 


Zoo 


ZccU 


rico, cm"^ 


Zoo 


•Zcell 


rieo, cm"^ 


Zoo 


5 


1782 


3.43 X 10^2 


933 


705 


1.09 X 10^^ 


381 


76 


9.27 X 10^^ 


44 


10 


1681 


4.02 X 10" 


1255 


710 


1.36 X 10i2 


536 


97 


1.21 X 10^^ 


75 


15 


1719 


1.22 X 10" 


1434 


744 


4.21 X 10" 


624 


109 


3.71 X 10^2 


93 


25 


1827 


2.79 X lO^o 


1654 


807 


9.87 X IQio 


733 


125 


8.48 X 10" 


115 


35 


1922 


1.07 X lO^" 


1797 


858 


3.82 X lO^o 


804 


137 


3.21 X 10" 


129 


50 


2038 


3.89 X 10^ 


1949 


918 


1.40 X lO^o 


879 


151 


1.14 X 10" 


145 


150 


2447 


1.73 X 10* 


2415 


1125 


6.36 X 10* 


1110 


180 


1.27 X 10^" 


178 



IV. DISCUSSION OF RESULTS 



Here, we have presented the results of theoretical stud- 
ies of the influence of the screened electric field on the 
charging of dust grains in a thermal plasma. We have 
applied the following approximations: (i) the cell model 
of dusty plasma for the description of the distribution 
of the electrostatic field; (ii) the quasichemical approach 
for the determination of the mean charge of dust grains, 
distribution of grains over charges, and variance of this 



distribution. 

The introduction of an electrically neutral cell contain- 
ing a grain allowed us to: (i) consider the circumstance 
that the emitted electrons remain in the vicinity of a dust 
grain, rather than move to infinity; (ii) describe satisfac- 
torily the distribution of the potential of the electrostatic 
field of a charged grain. Due to the introduction of a cell, 
we obtained a closed system of equations ((HI), ([II]), and 
(PS)) for the determination of the mean charge Z . 

The influence of the inhomogeneity of the distribution 
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Figure 5: 
2500 K, To 



Distribution functions pi{z) 
= 5 X 10~^ cm, and <; = 15 



and P2{z) at T — 
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Figure 6: Line corresponding to the values of rip and z (at 
T = 2500 K and rp = 5 x 10^^ cm), for wiiidi the condition 
g{z) = 0.9 is satisfied. Regions I and II correspond to g{z) < 
0.9 and g{z) > 0.9, respectively 



of the electrostatic field is significant for g{z) < 1. In 
Fig. [HI we demonstrate the line corresponding to g{z) — 
0.9 at the temperature T — 3000 K and the radius of 
grains rp 10^^ cm. Above this line (region I), we deal 
with densities rip and charges z, for which it is necessary 
to use the corrected formulas for the determination of the 
mean charge of dust grains. Below this line (region II), 
the influence of the screening effect can be neglected. 

We have shown that the screening causes an increase 
of the mean charge of dust grains as compared with the 
prediction of the theory omitting the screening effects. 
As the size of a cell increases, i.e., as the mean dis- 
tance between dust grains increases, result (|26p is trans- 
formed into the well-known Einbinder-Smith-Arshinov- 
Musin formula (f^ . 

For the grains of the radius rp = 10^^ cm, the differ- 
ences in the predictions of the two approaches are slight. 
This is explained by the fact that the grains of such sizes 
have very small charges ('^ 1), and the screening plays 
no significant role. 

In the experiments [s^l with grains Ce02 (Wq = 
2.75 eV), the mean density of grains Up varied in the 



limits (0.2-5.0) xlO^ cm~^, and the temperature T 
was changed in the interval (1700-2200) K. The mea- 
sured mean density of electrons ric was in the limits 
(2.5-7.2) xlO^" cm^"^, the mean radius of grains rp = 
4 x 10^^ cm, and the lower experimental bound of the 
mean charge Z w 500. 

For the grains with the indicated radius and the density 
7ip = 2 x 10^ cm^'^, the dimensionless radius of a cell 
<r sa 62. At the temperature T = 2000 K, the mean charge 
and the mean density of thermoemission electrons are, 
respectively, Z 442 and rie « 6.91 x 10^ cm~'^ according 
to ([^^ and ; the charge distribution variance is Zd ~ 
49 according to At the temperature T = 2200 K, 

we have determined Z w 560, ~ 8.76 x 10^ cm"-^, and 
Zd « 55. 

The calculated mean charge is close to the experimen- 
tal result {Z w 500). However, the calculated densities 
of electrons are less than experimental ones by one order. 
This fact can be apparently explained by the neglect of 
the influence of an ionized buffer gas in the calculations. 
Therefore, the comparison of the theoretical and experi- 
mental results is not quite proper. 

The study of the influence of the screening of a buffer 
gas will be performed separately. 

We note also that the essential excess of the mean 
charge Z over Zq is accompanied by the violation of the 
condition of applicability of the linearized Poisson equa- 
tion. The error of the calculated mean charge Z increases 
with Z , because, in this case, the contribution of the re- 
gion, where the nonlinear effects are of importance, in- 
creases as well. This will be studied in further work. 

The author thanks Prof. M.P. Malomuzh and 
Dr. V.I. Zasenko for the fruitful discussion of the results. 
The present work was partially supported by the Ministry 
of Education and Science, Youth and Sport of Ukraine 
(grant No. 0112U001739909). 



APPENDIX 
Electric Field of a Weakly Charged Spherical 
Grain in a Cell 

Small charges Z of a dust grain correspond to the condition 
^■0 ^ ^'c — f'p. In this case, the bulk density p of thermoemission 
electrons in a cell has a constant value p = — erie, where rio is 
the mean density of electrons in a cell, which is determined by the 
obvious realation 

^ _ 3 Z ^^^^ 
An rj? — rj^ 

The potential cp satisfies the Poisson equation 

Arff) = i-n-ene. (A2) 

In the dimensionless variables l[5}, Eq. (A2) supplemented by the 
boundary conditions (O has the solution 



i^if) = 1 + 



Z/Zo 



(A3) 



To within the designations, this formula coincides with that ob- 
tained in |16| in the case of the weak screening. 
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The strength and the energy of the electrostatic field of a weakly 
charged grain in a cell are given, respectively, by the relations 



1 Z^/Zq 

2 (<r3 - 1)2 



(A5) 



E{f) 



Z/Zo 



(A4) 



We note that, in the limiting case of small charges (without the 
screening), we have strictly -B(?) = 0. 
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